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In this study, a comprehensive view of a model crystal formation in a 
complex fluctuating medium is presented. The model incorporates Gaus- 
sian curvature effects at the crystal boundary as well as the possibility for 
superdiffusive motion near the crystal surface. A special emphasis is put on 
the finite-size effect of the building blocks (macroions, or the aggregates of 
macroions) constituting the crystal. From it an integrated static-dynamic 
picture of the crystal formation in terms of mesoscopic nonequilibrium ther- 
modynamics (MNET), and with inclusion of the physically sound effects 
mentioned, emerges. Its quantitative measure appears to be the overall dif- 
fusion function of the formation which contains both finite-size curvature- 
inducing effects as well as a time-dependent superdiffusive part. A quite 
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remarkable agreement with experiments, mostly those concerning investiga- 
tions of dynamic growth layer of (poly)crystaline aggregation, exemplified 
by non-Kossel crystals and biomolecular spherulites, has been achieved. 

PACS numbers: 05.20.Dd, 05.40.-a, 05.60.Cd, 05.70.Ln, 05.70.Np, 05.10.Gg, 
61.46.+W, 61.50.Ah 

1. Introduction 

The present study is devoted to modeling a formation of a soft (ordered) 
material known as the non-Kossel crystal ^j, as well as to propose how 
to model a polycrystalline aggregate termed the soft-matter spherulite 
The approach we are exploring throughout this paper assumes that the 
system of interest is a mesoscopic system P], by its nature combining both 
classical and quantum-mechanical properties, though we are offering here 
a description based on nonequilibrium statistical thermodynamics, taken 
suitably at a mesoscopic (molecular cluster) level jlj, i.e. a classical limit of 
the approach can preferentially be exploited [Hj. 

In general, examples of soft materials include biopolymers and charged 
polymer solutions, protein/membrane complexes, colloidal suspensions, gels, 
nucleic acids and their assemblies, to mention but a few. Most forms of con- 
densed matter - except of metals and ceramics (perhaps) jH] - are soft and 
these substances are composed of aggregates and macromolecules, with in- 
teractions that are too weak and complex to form crystals spontaneously; 
therefore, the task of how to grow crystals in complex environments still 
remains a real challenge T- 

Thus, the phrase 'soft condensed matter' has been coined, also for the 
clear reasons mentioned below. First, a striking feature is that slight per- 
turbations in temperature, pressure or concentration, or small (piconew- 
ton) external forces, can all be enough to essentially induce microstructural 
changes. Second, thermal fluctuations are almost by definition pronounced 
in soft materials, and in consequence, entropy is a privileged determinant 
of a soft microstructure formation, so that disorder, slow dynamics and ki- 
netics, and plastic deformation are the rules ^8,. Although soft materials 
have attracted engineers for ages, only recently have physicists taken an 
interest in such materials, and attempted to implement what is the essence 
of physics, that is to produce, at various levels of description, simple models 
that include the possible minimum information required to explain relevant 
features, cf. an example concerning proteins [U]. 

Nowadays, especially with the advent of single molecule spectroscopies 
are we managing to get quantitative data on things such as aggregation (and 
folding) effects in complex biomolecular environments, motion of molecular 
motors electron and proton transfer rates etc., that are sufficiently 
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reliable to formulate and test simple models. The availability of new exper- 
imental tools and simulation capability are urging physicists to apply their 
methods to structural biology, i.e. to allow the prediction of structure on 
the basis of known microscopic forces. 

A key issue is the fact that the environments in which a formation of the 
soft material happens are characteristic of pronounced spatial variations in 
dielectric constant, e.g. water and lipids jl2j . It is of major importance to 
realize that competition between interactions of different ranges results in 
different types of aggregation of molecules. The starting point should thus 
be about a discussion of the relative role of the various fundamental inter- 
actions in such systems (electrostatic, hydrophobic, conformational, steric, 
van der Waals, etc.), and what is their real impact on the aggregation of 
possibly ordered (crystal formation) and/or disordered types. The next fo- 
cus could be on how these competing interactions influence the form (and, 
topology) of soft and biological matter, like biopolymers and proteins, lead- 
ing to self-assembling systems of the type listed above 

In the underlying study, let us propose a way of modeling the soft matter 
objects named non-Kossel crystals jj], and in part also the biopolymeric 
spheroidal polycrystals, commonly termed spherulites. The way that we 
have chosen to achieve the above stated goal is based on the assumption that 
the biomolecular system, generally out of equilibrium, can be best described 
at a mesoscopic level, using the conception of mesoscopic nonequilibrium 
thermodynamics (MNET) [Tl], section 2. In sections 3 and 4, we will offer 
the thermodynamic-kinetic description of the soft material formation first 
by introducing its deterministic part (section 3), and second, by going to its 
stochastic, no doubt, more interesting part (section 4). In the deterministic 
part, we wish to place our emphasis on an unquestionable basis of each 
soft material formation, namely, how does the material formation depend 
upon the boundary effect, being typically curvature-dependent, and having 
also - in the case of spherulites - a nonequilibrium kinetic account readily 
involved in the overall process. In the stochastic part, in turn, we are 
lucky to support our mesoscopic view by a microscopic picture of a sub- 
process, essentially limiting the model material formation in the crystal 
boundary layer Here, we have in mind the physical fact that the 

growth of the objects that we model can thoroughly be controlled by the 
macroion velocity field of the crystal ambient phase nearby its interface with 
the growing crystal. Since it may lead to some very interesting physical 
consequences, e.g. a time-dependent viscosity effect ^B] (section 4), we are 
going to exploit this experimentally justified observation |17j in a sufficient 
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detail, looking at the macroion finite-size effect (FSE)^, and its impact on 
the overall crystallization process seen as a combination of static (curvature- 
oriented) and dynamic (constrained motion-oriented) effects in the crystal's 
growth layer ^5]. Concluding address, contained in section 5, closes the 
paper. 

2. MNET and its application to nucleation and growth 

phenomena 

Mesoscopic nonequilibrium thermodynamics (MNET) provides a general 
framework from which one can study the dynamics of systems defined at 
the mesoscale jl4(ll8j. and has been applied to analyze different irreversible 
processes taking place at those scales jl]. The formulation of the theory is 
based on the fact that a reduction of the observational time and length scales 
of a system usually entails an increase in the number of degrees of freedom 
which have not yet equilibrated and that therefore exert an influence on 
the overall dynamics of the system. Those degrees of freedom 7 (= {7i}) 
may for example represent the velocity of a particle, the orientation of a 
spin, the size of a macromolecule or any coordinate or order parameter 
whose values properly define the state of a mesoscopic system in a phase 
space. The characterization of the state of the system essentially relies on 
the knowledge of ^(7, t), the probability density of finding the system at the 
state 7 € (7,7 + d'j) at time t. One can then formulate the Gibbs entropy 
postulate ISl El in the form 



Here Seq is the entropy of the system when the degrees of freedom 7 are at 
equilibrium. If they are not, the contribution to the entropy arises from de- 
viations of the probability density P{^,t) from its equilibrium value Peqi'j) 
given by 



where AW = AW (7) is the minimum reversible work required to create 
that state is Boltzmann's constant, and T is the temperature of the 

^ In our approach FSE readily assigned to a macroion, will be emphasized twofold. 
First, it enters the nonequilibrium boundary condition of modified Gibbs-Thomson 
type (section 3). Second, it is involved in the change of viscosity near the crystal 
vs surroundings interface, the so-called size-dependent viscosity effect, cf. section 
4. Thus, the FSE is proposed to be a bridge between stochastic and deterministic 
parts of the approach. It will, for example enter the overall diffusion function of the 
complex crystal formation, D{R,t), see section 4. 




(1) 





GadomskiSiodmiak printed on February 2, 2008 



5 



heat bath. Variations of the minimum work for a thermodynamic system 
are expressed as 

AW = AE- TAS + pAV - fiAM + ctAS + . . . , (3) 

where (using a standard notation) extensive quantities refer to the system 
and intensive ones to the bath. The last term represents the work performed 
on the system to modify its surface S, whereas a stands for the surface 
tension. 

To obtain the dynamics of the mesoscopic degrees of freedom one first 
takes variations in Eq. 

5S = -kB I SP{j, t) In §^ d7. (4) 

focusing only on the nonequilibrated degrees of freedom. 

The probability density evolves in the 7— space according with the con- 
tinuity equation 

dP{l,t) _ dJ{j,t) 

dt ~ 97 ' 

where J(7, t) is an unknown probability current. To obtain its value, one 
proceeds to derive the expression of the entropy change, dS/dt, which follows 
from the continuity equation ^ and the Gibbs equation Q) . After a partial 
integration, one then arrives at 

-^ = -]Q-^Jsdi^a^, (6) 
where Js = J {"f , t)ln^^^^^ is the entropy flux, and 



is the entropy production which is expressed in terms of currents and con- 
jugated thermodynamic forces deflned in the space of mesoscopic variables. 
We will now assume a linear dependency between fluxes and forces and 
establish a linear relationship between them 

Jh.,) = -feLwAL^V (8) 



where L = L{'y) is an Onsager coefficient, which depends on the mesoscopic 
coordinates 7; in general, it also depends on the state variable P{^,t). To 
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derive this expression, locality in 7— space has also been taken into account, 
for which only fluxes and forces with the same value of 7 become coupled. 

The resulting kinetic equation then follows by substituting Eq. (jSJ back 
into the continuity equation 

= A (d(. t)Ph t)^ln^] (9) 

where we have defined the diffusion coefficient as D{^,t) = This 
equation, which in view of Eq. @ can also be written as 

mi.t) _ a /oh,t)^ + :^^Ph.t,) , (10) 



dt V d'y ksT O7 

is the Fokker-Planck equation accounting for the evolution of the probability 
density in 7-space. 

2.1. Nucleation processes 

The expression of the nucleation rates can be obtained from the proposed 
formalism 5 . To this purpose one has to interpret the nucleation (or in 
general any activated process) as a diffusion process along the mesoscopic 
coordinate describing the state of the system, the embryos, at short time 
scales in between the metastable and the crystal phases. The reaction rate 
or diffusion current can be written in terms of the fugacity z = exp/z/fesT 
as 

1 dz 

J{^,t) = -kBL- — . (11) 
z (77 

The current can also be expressed as 

dz 

J{j,t) = -D{^,t) — , (12) 

where D{^,t) = ksL/z represents the diffusion coefficient. As a first ap- 
proximation we assume that D('y,t) = D = const, and integrate from the 
initial to the final position, obtaining 

J ^ r jdj = -D{Z2 - z^) = -D{exp - exp -^). (13) 
J^i I^bT KbT 

This equation can alternatively be expressed as 

J = Jo {1 - e^/'-^) , (14) 
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where J is the integrated rate, Jq = Dexp{fii/kBT) and ^ = //i — ^2 is the 
affinity. We have then shown that a mesoscopic thermodynamic analysis 
may lead to the formulation of the nonlinear kinetic laws governing nucle- 
ation processes. Nucleation processes in which the embryos are embedded in 
an inhomogeneous bath can also be studied by means of the theory proposed 
0. 



2.2. Agglomeration and growth processes 

MNET can in general be used to study kinetic processes taking place 
in mesostructures such as the nanostructure arrays The size of these 
structures is in between those of single particles and macroscopic objects. 
They carry out assembling (clustering) and impingement (pattern forma- 
tion) processes. They may also diffuse in a thermal bath, be convected by 
external flow and be affected by external driving forces as those acting in 
extrusion, shearing, injection processes involved in the mechanical process- 
ing of a melt. The growth rate of the agglomerates can be determined by 
the formalism previously presented in which the mesoscopic variable is the 
volume of the grain. One formulates the Gibbs entropy postulate in terms 
of the probability density P{v,t), where now the degree of freedom 7 is 
V - the volume of the grain or molecular cluster. Proceeding as indicated 
previously one obtains the growth rate 



TP{v,t) 



(15) 



Interpreting $ as an entropic potential 22 and assuming that the volume- 
dependent Onsager coefficient L{v) follows a power law of the type v'^, where 
6 = 1 — 1/d, with d the dimension of the system, one obtains the expression 
of the rate [El 1221 

J(v, t) = -aav'-'P{v, t) - Dav'^^^ (16) 

where aa and Da are reference constants ^Sl- This expression constitutes 
the Louat-Mulheran-Harding {L — M — H) law I2S1 proposed heuristically 
by assuming that the drift is due to surface tension effects, considering only 
the mean curvature. The procedure can be generalized to the case in which 
the Gaussian curvature is important which occurs at small sizes of the grains 

1221. 

In the offered modeling we would explore another degree of freedom 
which is i?, the radius of a molecular cluster or a crystal. It appears natu- 
rally in the modeling being of the form of Fokker-Planck type ^ . 
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2.3. An example: 2D modeling of a crystal growth in complex environment 

and its MNET features 

A multigrain growth has also been considered as a discretized process 
in terms of MNET in both space and time domains jH]. The model is meso- 
scopic in the sense that each cell can contain (i) a liquid unit of phase I, (ii) 
a crystal unit, (iii) a small particle embedded in some melt or (iv) a solid 
large particle. Initially, each cell of the lattice contains either a large parti- 
cle with some probability P;, or a small particle embedded in some liquid I 
with a probability Rj, or a liquid unit with a probability 1 — Pi — Pg. Nucle- 
ation is induced by simultaneously turning a number n of cells into initial 
solid units, supposed to be randomly dispersed in the initial melt. Since 
the process occurs over long time scales, thermodynamic (i.e. equilibrium- 
like) quantities can be mapped into multigrain growth probability rules, as 
follows. At each growth step, all i cells containing some liquid phase, i.e. 
/-cells and cells with a small particle, in contact with mesoscopic cells are 
selected. The probability Pi to grow the phase on the cell i is given by a 
classical thermodynamic argument as Pi ~ exp(—AGi/kBT), where AGi is 
the gain of free energy. Usually, it can be decomposed into two terms: a bulk 
contribution depending on the driving force and a local surface contribution 
which is proportional to a chemical bond energy E^, cf. Eq, 

The algorithm on which the discrete modeling is based, combines a mech- 
anism for pushing and/or trapping of particles, a chemical reaction, and 
crystal growth kinetics. The model is an adapted Eden model Pl]. The 
change in energy enters in the Boltzmann factor, cf. Eq. 0. It is con- 
sidered to be anisotropic and depending on the number of neighbors. This 
leads to rugged surfaces and facets. Moreover, interaction with unreactive, 
or reactive particles has been considered as in the theory of particle trap- 
ping or displacement along growing interfaces, elaborated by Uhlmann et 
al., and called the U — C — J model j^Sj. The particles are trapped or pushed 
by the mobile solid/liquid interface depending on (i) the surface tensions at 
solid/liquid and particle/liquid interfaces, (ii) the particle (impurity) size 
and (iii) the growth velocity of the solidifying front. (A similarity between 
the U — C — J and what we have proposed in sec. 4 has to be noticed, cf. 
Fig. 3 therein.) For a fixed particle size (see, the FSE effect involved in the 
modeling presented in sections 3 and 4), the particles are supposed to be 
trapped by the front if the growth velocity of the interface is higher than 
a critical value, a physical situation that suits perfectly our type of model- 
ing. This critical value is controlled by the particle size and the interfacial 
tension energies j25j . 

Since, while growing the phase in an entropic milieu, the bulk contribu- 
tion is roughly constant in the system for isothermal conditions, only the 
local surface contribution is needed for measuring the probability of growth, 
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i.e. Pi ~ exp{—gnn^i) where Ni is the number of nearest-neighboring (nn) 
units belonging to the same grain and gnn = Ef, x nn/kBT. This expresses 
an anisotropic locally preferred kinetic growth along the main directions. 
For high positive values of square- like grains are growing j^H]- For 
large gnn, smooth grain faces are obtained. 

^From a kinetic point of view, it is accepted that the grains grow anisotrop- 
ically with a four- fold symmetry along the a — b planes, the edges of the 
grains being oriented along the (10) and (01) crystallographic directions. 
This has been implemented in the choice of the parameter value gnn- 

In order to simulate best the chemical reaction process, the mesoscopic 
cells for which the nearest and next-nearest neighbors do not contain any 
particle are excluded from the growth cell selection. Indeed, the growth of 
the phase on these mesoscopic cells is assumed to be improbable because of 
the local deficiency of the basic component around these cells. Taking into 
account the probabilities Pi of all possible growing cells, one specific growth 
site is randomly selected. 

The simulations [23 UHl IB] put into evidence the effect of the particle (no 
matter whether a 'domestic' or testing particle, or some impurity particle) 
size distribution itself on the resulting microstructure, which is of special 
interest of the underlying study. Even though, the MNET-based model 
considers only two different particle (impurity) sizes, the model predicts that 
the refining of the impurity particles leads to better samples. It should also 
be underlined that a discrete MNET-oriented modeling can be thought of 
as a useful working extension of many of its continuous versions jl4l 1181 122j . 

Some extension from 2D to 3D is rather straightforward but not too 
trivial: Just as difficult as going in any simulation of Eden-type systems 
when counting sites during spreading must be done with care, see [21], and 
refs. therein. 

3. Deterministic part of the approach to a sphere growth 
controlled by mass-convective fluctuations: the role of 
curvatures and nonequilibrium boundary effects 

In this section, we consider two physically relevant cases of the bound- 
ary condition (BC) prescribed for the growing, e.g. protein (non-Kossel ^), 
spherical crystal. In the first case, we propose an essential modification of 
the curvatures' contribution to the BC. This includes the Gaussian curva- 
ture. This condition is of equilibrium type. In the second case, we propose 
" Goldenfeld-type" condition j27j which is of nonequilibrium type. In the 
equilibrium case we would like to offer a more detailed analysis, whereas 



^ Non-Kossel crystals are defined as complex structures with several molecules per unit 
cell in inequivalent positions 7 . 
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in the nonequilibrium case, because the analytic solution to the problem 
cannot be found, we only discuss the asymptotic limit in a brief way. 

Let us assume that initially at t = the growing object is an ideal sphere 
of radius Rq and the density of the sphere is C. At time t > the radius 
of the growing sphere is equal to -R = R{t). Prom the mass conservation 
law ^2] an evolution equation that gives the speed of the spheroidal crystal 
formation arises, and looks as follows 

dR Csjl/R) 

where Cg = Cs{l/R) is the curvature-dependent surface concentration, and 
v{R) is the velocity of incoming macroions, both taken at distance R from 
the sphere center. The concentration c^, prescribed at the boundary, is 
derived under the assumption of local thermodynamic equilibrium at the 
boundary. It has the form of the well-known linearized Gibbs-Thomson 
relation 

c, = co(l + 2r4), (18) 

provided that T is the so-called Gibbs-Thomson or capillary constant, which 
is usually of the order of Wnm for lysozyme crystal 0, cq is an equilibrium 
concentration for the planar surface, practically for R^ Rq. Based on the 
knowledge of formation of droplets viz (vapor) condensation processes, one 
typically defines T as j^H] 

where m is the mass of the vapor-phase atom (unit) and pd stands for the 
density of the droplet's material, being here identified with C of Eq. (|17|) : 
a and T have their usual meaning (see above). The above definition can 
also be adopted to our case fT5]. 

The solution of Eq. H17() with Gibbs-Thomson boundary condition, Eq. 
Pll. for v{R) 

— '^mi — const, reads jl2j 

R -\- 2r 

R-Ro-{Rc + 2T)ln = aoVmit, (20) 

Kq + zi 

and its large time asymptotic solution becomes 

Rr^t. (21) 

Here 

^0 = (22) 

O — Co 
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Fig. 1. Schematic picture of the mean (global) and Gaussian (local) curvature 
problem. Twice the mean curvature of the growing sphere of radius R equals 
and is a global curvature measure. Upon magnifying the visualisation box (VB) 
a local picture appears, at a given time instant t and packing conditions "q (see 
caption to Figure |2J), for which the Gaussian curvature of the building block viz 
the spherical molecule of radius r,„i reads The overall curvature contribution 

to the equilibrium BC of Gibbs-Thomson type, Eq. H18|l . reads ;r^-^, - Tolman 
parameter |23l, and leads to a modified BC, Eq. H18I) . of importance for crystal 
composed of finite-size (or non-point like) building blocks. 



and 

Rc = 2fTor (23) 

is a critical nucleus' radius, fio is an equivalent of the bulk supersaturation'^. 

As the first essential modification (equilibrium type) of the curvature 
contribution to the BC we propose to include the Gaussian curvature. 
Schematic picture of the mean (global) and Gaussian (local) curvature prob- 
lem is presented in Figure ^ This correction can be important, because of 
large linear dimension of the growth units in comparison with linear di- 
mension of the other solution components, so the local curvature can play 
important role, especially when we deal with macromolecular islands built 

^ Except that the original Gibbs-Thomson condition involves a firm curvature- 
dependent contribution, see Eq. II18II . it is rather hardly applicable to crystallization 
under high bulk supersaturation. This fact imposes some limits to the magnitude of 
the overall chemical-potential based driving force of crystal growth as a whole |7ll31|. 
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up of a few units. Now the BC takes the form 

cs = co(l + 2ri + ^i) = co(l + 2fi), (24) 

tX It 

where 

r = r + -^ (25) 

^' mi 

and l/{rmiR) is a Gaussian curvature, Vmi is the radius of the crystal build- 
ing unit (for lysozyme protein r^i ~ 1.5nm) and 6t is a Tolman length. 
(5r is defined as the difference between the radius of the surface of ten- 
sion and the radius of the equimolar dividing surface (for lysozyme protein 
St ~ 3.5nm) and depends on the packing coefficient of the crystal. The 
Tolman length can be related to the superficial density at the surface of 
tension [^. Note that by postulating Cg in the form of Eq. H24|). with 
6t involved, we somehow induce, for a given time instant t, elastic effects 
at the crystal boundary. This is because 6t is also a measure of rigidity 
for bending of a curved piece of the crystal boundary, a physical fact well 
known to those studying biomembrane formation j^. Thus, we may, also 
by means of having 6t involved in the BC, arrive at elastic contribution 
to the crystal growth This situation resembles a well-known fact that 
an improperly placed boundary molecule, or some impurity, exerts an ad- 
ditional strain-stress field in its very vicinity j31j . 

The solution of Eq. with BC prescribed by Eq. (PH), with is of 
the same type as the one given by Eq. (|^n)) but now 

Rc = Rc = 2aof (26) 

has to be inserted, cf. Eq. if^ . 

The dependence of R{t) is presented in FigureEJ The lower-right inset to 
Figure 121 shows differences in behavior of the surface concentration with and 
without Gaussian curvature correction to the surface concentration, while 
the upper-left inset points to early-stage differences in tempo of crystal 
formation expressed in terms oi Rvs t dependence. These differences stand 
also over longer time period but they cannot be seen as separated curves for 
time interval (1.5 hour) chosen to fit our model to experimental data [6'2\ . 
We can see that Gaussian curvature influences (speeds up) the growth rate 
just by increasing the local surface concentration, especially for the early 
stage of the growth process when R/rmi is relatively small. 

Our results adjust well to the experiment. Chow et al. |32j observed 
constant growth rate of the lysozyme spherulites, what is characteristic 
when the growth is controlled by surface kinetics rather than by volume 
diffusion. In the experiment the crystal density is 25 times bigger than 
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Fig. 2. The time dependence of the crystal radius R{t) for original H18|l and modified 
(|24|) BCs of Gibbs-Thomson (and Tolman) types taken for three different crystal 
structures (cubic close packing rji, hexagonal lattice 772 and cubic lattice 773). The 
lower-right inset shows differences in behavior of the surface concentration with and 
without Gaussian curvature (finite-size or Tolman type) correction to the surface 
concentration, while the upper-left inset points to early-stage differences in tempo 
of crystal formation. These differences stand also over longer time period but 
they cannot be seen as separated curves for time interval chosen to fit our model 
to experimental data (C = 500 [arbitrary units], cq — 1 [arbitrary units], T = 



1 • lO-^TO, Rc = 4:- IQ-"™, i?o = 5 • lO-^T 



= 2.3 • 10"^m/s). 



lysozyme concentration in the solution. In our calculation supersaturation 
is equal 500. This discrepancy implies that the narrow interfacial region 
(double layer) designed for our model is definitely a protein-poor region. 

In the second case we propose "Goldenf eld- type" correction [23- Now 
we assume that the surface is away from local equilibrium and deviation 
from this state is proportional to the growth velocity of the interface 

c, = co(l + 2f4-/?G^), (27) 
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where (3g is a positive kinetic coefficient. Solving Eq. (|17j) with nonequi- 
hbrium boundary condition, Eq. (|27() . one gets for long times the linear 
solution of exactly the same type than that given by H21|) which can re- 
flect the asymptotic temporal behavior of spherulites. It can be shown that 
when increases the growth becomes slower. At early stage of the growth 
process the influence of (3g on R{t) is strongly nonlinear. For long times, 
the growth is linear in time and Pg influences only the growth rate (slope 
of the curve, cf. Figure [21) 2 . 

Looking at the Eq. (|17jl we can see that the growth process depends 
not only on the supersaturation of the solution. It turns out to be very 
important that the macroion velocity, which depends on the electrostatic 
interactions of the macroion with the solution (we have to mention that the 
solution is an electrolyte), also depends on the viscosity of the solution, and 
moreover, on the size of the growth units. It is sometimes possible that the 
growth unit will be made of few macromolecules, so that in the viscoelastic 
complex environment that we actually describe, the aggregates (made up 
of 2 — 5 macromolecules) will have different velocities, especially when com- 
pared with that of the single macroion because the diffusion coefficient can 
now be aggregate-size dependent |33| . This experimentally justifled obser- 
vation 34. has its profound impact on the so-called intrinsic properties of 
the viscoelastic solution ^Bj, making them time-sensitive, thus imposing an 
intrinsic time scale, presumably powerly "correlated" with the observation 
time scale jSSl. It will be demonstrated in the subsequent section. 



4. Sphere growth controlled by convective fluctuations: the role 
of time-dependent viscosity and flnite-size effects 

In accordance with the considerations of the previous section, here we 
will complement the deterministic description of the kinetics of crystal 
growth by incorporating thermal fluctuations of the bath through fluctu- 
ations of the velocity of the incoming macroions (or growth units) j36j . 
and relating them with the finite size of the particles at the locally highly 
concentrated regime. In this regime, the viscoelastic properties of the so- 
lution become important and must also be incorporated into a consistent 
description. 

To perform this description, we will adopt an effective medium theory 
j37j in which the diffusion of a "test" macroion through a saturated solu- 
tion can be analyzed in terms of the single-particle distribution function 
f{f,u,t), depending on the position r and instantaneous velocity u of the 
macroion. In accordance with mesoscopic nonequilibrium thermodynamics 
formalism, in this case the entropy per unit volume, s{t), can be expressed 
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in terms of the Gibbs entropy postulate |SZ1 OHj 

s{t) = si^-kB I fln—du, (28) 

where si^ is the entropy at local equilibrium, cf. Eq. and we have 

introduced the fugacity z = a f{f,u,t), with a the activity coefficient char- 
acterizing the interactions of the macroion. zie is the fugacity at local equi- 
librium. The activity coefficient a will contain two contributions since the 
test macroion interacts with both, the solution of macroions (mostly, in 
the vicinity of the crystal boundary) and with the crystal Oc, see Figure 01 
The activity of the solution can be expressed as = e'^ p/ksT ^ with c~^p 
the compressibility factor in which p is the excess of osmotic pressure |37) 
and the concentration field is defined by c(r, t) = f f{f, u, t)du. The explicit 
expression of c~^p can be obtained in the context of different theories on 
electrolyte solutions jSH]- The crystal activity Oc can, in general, contain ef- 
fects related with entropic, energetic or geometrical constrains, as indicated 
in section 2, 11]. 

Following the procedure of MNET already outlined, one may show that 
the Fokker-Planck equation describing the evolution of f{f,u,t) is [371 l4Uj 



^ + V . (u/) - ksTV [In a, + In a,] . 1^ = ^ . /3(t) 



\ kBTdf 
uf + 



m 



du 



(29) 



where the third term on the left-hand side of the equation contains the 
forces acting on the macroion. The right-hand side of this equation con- 
tains the contributions of Brownian motion of the macroion in which the 
time-dependent friction coefficient [3{t) introduces memory effects j411 142j . 
As mentioned in the previous section, these effects are related to the vis- 
coelastic properties of the solution (time-dependent viscosity), and modify 
the intrinsic time scale of the particle dynamics. 

Several models for the friction have been used in the literature of model 
crystal growth j^H]- However, since we are interested in analyzing the effects 
of the finite size and inertia of the macroions, then we must take into account 
that in this case a general expression for the friction coefficient is j4U[ I37j : 
I3~^{t) oc [1 + {t£)UjY''']~ , with the symbol denoting the inverse 
Laplace transform, uj the frequency and td = -ff^ a characteristic diffusion 
time, with Dn = ^ — the diffusion coefficient of the macroion at infinite 
dilution, r^i is the macroion radius and r/j^ the viscosity of the solvent 
at infinite dilution. In this relation, the exponent 7c may in general be a 
function of the macroion concentration at the bulk . In order to make 
an analytical progress, for simplicity one may use the expression of /3(t) 
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© CI ion Na* ion 



Fig. 3. Naive picture of the vicinity of the crystal surface. Building unit (posi- 
tively charged lysozyme protein), which is surrounded by neutral (water dipoles) 
and oppositely charged particle (counterions) , performs random walk (RW) with 
a constant velocity in an electrostatic double layer. We can see that the crystal 
surface is not smooth and the local (Gaussian) curvatures are shown (see caption 
to Figure^. In the left top corner a triplet has already been formed. 



obtained in the context of the generahzed Faxen theorem [401 143j 



exp I — I Erfc 
\td) 




(30) 



where /?o = Q'n'r],id'>'mi is the Stokes friction coefficient. Notice that Eq. (|5n)) 



and that 



is the inverse Laplace transform of /3 ^{uj) = Pq^ 1 + (tdw) 
at times t << td, it reduces to I^q^ \J^f-- 

To obtain the velocity v{t) of the macroions entering into the growth rule 
H17|l in which, however, the macroion velocity has been assumed constant, 
see legend to Figure |^ one may use Eq. (|29|) to obtain the hierarchy of 
evolution equations for the moments of /, and from it construct the equation 
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(if* V 



(31) 



which describes the evolution of the macroion in the diffusion regime. Here, 
we have defined the collective diffusion coefficient = kBT[3~^ {l + ■ 
Particular expressions for this coefficient can be modelled within the scope 
of, for example, Kirwood theory of solutions. In a mean field approxi- 
mation, the pressure p can be approximated by an expansion of the form 
pksT [c — B{T)(P'\ , where B{T) is a virial-type coefficient incorporating the 
specificities of the pairwise interactions among macroions. In general, B{T) 
may depend on the radial distribution function, and then on spatial corre- 
lations. Substituting this expansion into as = pI^bT ^ taking the loga- 
rithm and the corresponding derivative of the definition of Do one obtains 
the time-dependent diffusion coefficient Dc = ksT {1 — B{T)) (3~^{t). li is 
interesting to notice that T(l — B{T)), resembling an "effective" tempera- 
ture j37j . implies that macroion diffusion has a nonlinear dependence on T, 
which could be relevant in the crystal structure, as was already mentioned. 
From Eq. H31() . it follows that in the diffusion regime the velocity of the 
macroions will satisfy a Langevin equation of the form 



dr{t) 
dt 



v{t)=(3-\t) \Fd{t) + F\t) 



(32) 



where for simplicity we have considered the unidimensional case and F^{t) 
constitutes a random force due to thermal fiuctuations. Moreover, we have 

c{r,t)dr, due to the 



B 



rpdlnoc, 
dr 



defined the deterministic force F(i{t) = J 
interaction of the macroion with the crystal. Notice that the fiuctuating part 
of the velocity satisfies the usual conditions of zero mean and {v{t)v{ti)) = 
p-\t-ti) _41.. 

In a first approximation, one may assume that in the vicinity (double 
layer) of the crystal surface, the electrostatic force (attracting macroions) 
is constant, = Fq, ^21- This assumption implies that Uc depends upon 
the minimum reversible work necessary to move the macroion a certain 
distance Ar, Uc = e^o^r/ksT ^ ^ taking into account the 

above assumptions, near to the crystal, the growth rule ((T7jl can be written 
in the approximated form 



dR{t) 
dt 




Fo + F«(t) 



(33) 



where = . In the present approximation, an can, for instance, be 
expressed in the form 
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^R = '^o J-, (34) 

^ iJ 

after using Eqs. (PU, (EHl) and Note that for t >> 1 ctr ^ o-q holds. 

Moreover, notice that for Cs we have to take here, depending on what we wish 
to model, either Cs from H24|) (equilibrium BC: non-Kossel crystals) or the 
one from ()27() (nonequilibrium BC: biopolymeric spherulites) . Especially, 
the latter seems to be a challenge for extensive numerical modeling, which 
is, however, left for a future task. In both cases mentioned, however, the 
FSE account is thoroughly manifested, cf. Figs 1-4. 

Now, in order to show the relation existing between the present approach 
with that of Sec. 2, for simplicity we will consider the case in which the 
deterministic part Fq can be neglected. As a consequence, Eq. becomes 
a stochastic equation with multiplicative noise which is equivalent to the 
Fokker-Planck equation 



dP{R,t) _ d 
di ~ 'dR 



. ^ dR keT dR ^ ' 



(35) 



where P{R, t) is the probability density in the i?-space with R standing for 
the radius of the spherical crystal. Moreover, we have defined^ the potential 
<I> = In ctr and the diffusion coefficient 

D{R,t) = DR{t)al, (36) 

where 

DR{t) = dJ— (37) 

is the time-dependent diffusion coefficient of the incoming macroions related 
with I3~^ through jHEl 

Djiit) = ksT f' (3-\h)dh. (38) 
Jo 

Figure 01 shows the time-dependent part of the reduced diffusion coefficient 
as a function of the reduced time. The solid line represents the integral 
of Eq. (jSOJ, whereas the short-dashed line is represented by Eq. (|n7|) . see 

* In another paper Q we have defined "I> = ksTlnaR, i.e. in a Boltzmann-like form. 
It automatically implies that the time-dependent part of Eq. I).36|l . DR{t), has to be 
defined without ksT as a prefactor absorbed in Dq. It, however, destroys an Stokes- 
Einstein type form of Do, which is not the case of the present paper in which this 
form is used explicitely, cf. caption to Figure |^ 
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Eq. ■ The long-dashed hne represents a behavior for comparison. 
From the solid line it follows that at short times the crystal grows with 
a super diffusive behavior whereas at larger times, Dn tends to a constant 
value, implying that the crystal radius tends to a maximum value, i.e. to 
a final equilibrium state. It coincides well with both the deterministic view 
of the process, and above all, with its experimental realizations |151 131j . 
It is interesting to realize that formally both the growth speed and the 
constrained motion of a 'representative' macroion, feeding the crystal go 
in a superdiffusive way. Thus, it seems that both sub-processes manifest a 
kind of synchronous temporal behavior. 

Finally, notice that from equation ()35() it follows that the deterministic 
part of the current is 

Comparing this expression with Eq. (|12|) , one may conclude that (Jr plays 
the role of a fugacity in i2-space whereas ^> the role of a chemical poten- 
tial containing geometric restrictions due to the boundary conditions, as 
discussed in the previous section ^Hj. 



5. Summary and conclusions 

In this paper, we have proposed an integrated static-dynamic picture 
of model (spheroidal) crystal growth in a complex milieu, preferentially 
of electrolytic type. The complex milieu we have in mind here could, in 
general, be a two-component system, containing some impurities and/or 
(charged) additives. The spectrum of examples can range from metallic 
polycrystals (L — M — H model, sec. 2.1 and 2.2) via superconductors^ 
(U — C — J, sec. 2.3) until the soft-condensed matter objects that are: 
non-Kossel crystals and (bio)polymeric spherulites. In fact, the latter suits 
best our type of modeling since it really manifests the FSE property due 
to pronounced sizes of the macromolecules constituting the (dis)ordered 
macromolecular cluster. The goal of readily incorporating the FSE prop- 
erty into the growth rules governing spheroidal crystal formation has been 
successfully achieved at both static as well as dynamic level of description 
- therefore our integrated model is called static-and-dynamic. At the static 
level of description the goal has been achieved by proposing a modification 
of the well-known Gibbs-Thomson (equilibrium) formula or by additionally 
including its nonequilibrium modification that we attributed to Goldenfeld, 

^ YBCO-\ike superconducting ceramics grown near a peritectic point i.e., when a solid 
is in contact with a liquid and there is incomplete reaction |25| 
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for 



Reduced diffusion coefficient -rr as a function of the reduced time — 
Dq ~ 10~^CTO^s~^ and td — 10~®s. The solid line represents the integral of Eq. 
H3U|I , whereas the short-dashed line is represented by Eq. 1)3 7|) . In the deterministic 
case, the integration of Eq. (|33l) can be compared with that of Eq. H20|l with the 
parameters F and i?c, from which follows that at sufficiently large times the left- 
hand side behaves as R, whereas the right hand side will be proportional to Eq. 
(|38|) . Since the solid line tends to a constant asymptotic value, then from the 
figure it follows that for sufficiently long times, the radius of the crystal tends to a 
maximum value. In contrast, (3^^ implies a "pure" superdifFusive growth. The 
long-dashed line, corresponding to a behavior, has been included for comparison. 
It can be obtained form fi~^(t) oc L^^ [1 + (r^jw)''''] ^, when the exponent 7c = 
1/4. 



who likely did it for the first time almost twenty years ago but for a diffu- 
sive field feeding the growing object |27j. In this study, however, we have 
concentrated on the equilibrium BC with the Tolman-type correction pro- 
posed, emphasizing a pivotal role of the Gaussian curvature in the realistic 
view of the process that we have offered. Therefore, the nonequlibrium type 
BC is dealt with in a rather sketchy or qualitative way, but nevertheless, 
such a possible modification is worth mentioning, at least for some inter- 
esting future task. Some undoubtful integrity of the picture, or its quite 
comprehensive character, at which we have finally arrived, can by no means 
be accomplished when no complementary stochastic part had a chance to 
appear. This is, in our opinion, the most worth-emphasizing part of the 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



R[m] X10-' 

Fig. 5. Mean course of the overall diffusion- type behavior of the non-Kossel (model) 
crystal growth. The curvature effect and FSE are largely pronounced for Ic. ~ \ 
and some consecutive time steps t. The same type of temporal behavior has been 
proposed in |36j for algebraic correlations with exponent 7c = ^. Solid lines are 
plotted on the basis of Eq. H3fci|) . whereas the dashed lines on its equivalent with 
Rc = Rc and r = r (for details see caption to Figure 0J. 



description, very responsible for unveiling its dynamic aspects. It could also 
be worth-presenting when forseeing a necessity to go beyond some analytic 
modeling El EH] > and to propose thoroughly computer simulations, 

carried out in a complementary way [SH^HHni- Thus, the integrated static- 
dynamic picture of the process is offered in sec. 4. It tells us that: (1) the 
role of curvature(s) is important; (2) the role of FSE is equally important, 
and anticipates the approach to be quite realistic; (3) the role of velocity 
fluctuations in time appears to be crucial in setting up properly the picture: 
a basic result of sec. 4 points to the constrained motion of the macroion 
to be superdiffusive, with the exponent 1/2 for an unidimensional approxi- 
mation offered therein. Thus, and in particular, a specific goal to propose 
a firm mathematical form for algebraic correlations, somehow underscored 
in another study j36j has been achieved too. The account of correlations in 
space is not explicitely taken into account in this study. Implicitely, how- 
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ever, it is involved in it by proposing Eq. with the collective measure, 
Dc, from which further a Langevin-type as well as Fokker-Planck type de- 
scriptions for model crystal growth (within the integrated picture) arise. 
What above all comes out from the proposed integrated picture, however, 
is a clear dynamic measure of this integrity, represented by D{R, t) of Eq. 
H36|) . It is seen in a picturesque way in Fig. 5 which tells us that the pro- 
cess goes naturally in the way that for its final (mature) growing stages 
the overall diffusivity of the system, due to the above mentioned physical 
factors (FSE, curvature(s), fluctuations) slows down, this way likely pro- 
moting order against disorder in the complex entropic milieu from which 
the object emerges. Also, a difference between Gaussian and non-Gaussian 
curvature effects on the diffusive behavior of the system clearly pronounce 
in the course of time. Finally, let us underline a remarkable feasibility of 
the presented type of modeling to serve for fitting to some experimental 
curves, or to relate the proposed theory to the experiment. An example of 
such an ability has been offered by Fig. 2, i.e. the case of biopolymeric 
spherulites, for the kinetic coefficient, 0, i.e. when the nonequilib- 

rium BC gets finally equilibrated. Other fits can be done for non-Kossel 
crystals as well iSljilSj. Last but not least, let us juxtapose some features 
of MNET, because this thermodynamic-kinetic formalism enabled to offer 
such an integrated view of the process. These are [TH 151 IIHIH^ : 

• MNET provides a description of kinetic proceses taking place at the 
mesoscale (and involving mesostructures) in which fluctuations are 
important. 

• The theory provides in general expressions for the rates (nucleation, 
growth, reaction, etc.) which are in general nonlinear functions of the 
driving forces. They are derived in very general conditions even when 
mesostructures are embedded in an inhomogeneous bath. 

• The theory provides kinetic equations of the Fokker-Planck type for 
the evolution of the probability distribution. From them one can de- 
rive expressions for the moments of the distribution which are related 
to physical quantities which can be compared with experiments. 
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